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I. INTRODUCTION 

Quantum measurement is usually taught in textbooks 
as an instantaneous process. However, in nature, all pro- 
cesses take a finite time. Physical projective measure- 
ment (or wave-function collapse) must therefore happen 
over some time period, and is often a sequence of weak 
measurements, run for a sufficiently long time. Weak 
quantum measurements are characterized by an intrinsic 
uncertainty about the state of the measured system. In 
the parlance of detector physics, this is equivalent to the 
statement that the signal cannot be confidently distin- 
guished from the noise without a sufficiently long integra- 
tion time. At some intermediate time, the eigenstates of 
the measurement operator can only be assigned a value 
with some confidence, determined with the probability 
of a given realization of the detector output (similarly to 
classical Bayesian inference). 

In solid state systems, the typically weak coupling con- 
stant between the quantum system and measuring appa- 
ratus imply that weak measurements with long measure- 
ment times are the norm. While this is often frustrating 
to experimentalists who wish to preform projective mea- 
surement for quantum computation purposes, we view 
this situation as an opportunity to discover and imple- 
ment ideas in quantum measurement that are qualita- 
tively different from standard projective measurement. 
Although many results in this paper are abstract and ap- 
ply to many different physical systems and detectors, the 
results will be discussed in terms of solid state physics. 
Quantum detection in mesoscopic structures began with 
the "controlled dephasing" experiments of Ref. [l| and re- 
lated theoretical works which has continued to be an 
active area of research)^ The particular mesoscopic struc- 
ture we shall focus on presently is a quantum point con- 
tact (QPC) detector measuring a double quantum dot 
charge qubit (DD), a system that has been extensively 
investigated, both theoretically^ and experimentallyi^ 

A generic problem that arises if one is interested in 



making a projective measurement, made out of many 
weak measurements, is that the dynamics from the 
Hamiltonian evolution combines in a nontrivial way 
with the measurement dynamics. The way around 
this problem is with quantum nondemolition (QND) 
measurements The QND scheme employed in this pa- 
per is that of kicked QND measurements ''■^^^MSjiLi^ on a 
qubit. This idea was discussed for the QPC in Ref.0, and 
was inspired by a similar rotating QND scheme of Ref. 13. 
By turning the detector on and off in a time scale much 
faster than the Rabi oscillation period, the detector gives 
a little information about the quantum state. The qubit 
is then allowed to make a full Rabi oscillation before the 
next weak measurement, so the measured observable is 
static in time from the perspective of the measurement 
device, effectively turning off Hamiltonian evolution. In 
this fashion, the information contained in the qubit is 
teased out over many measurements, or detector kicks. 
Kicked QND measurements have many advantages that 
recommend them as a technique of choice both for theo- 
retical treatment, as well as for experimental implemen- 
tation: 

• Theoretically, the kicked measurements may be de- 
scribed with a non-unitary quantum map. The map 
is discrete in the time index, but the measurements 
are weak, not projective. The quantum map for- 
malism allows for a technically simple treatment 
of the combined Hamiltonian and measurement- 
induced dynamics in a global manner, in contrast to 
continuous measurement analysis using conditional 
differential Langevin equations. 

• The kicking mechanism is convenient, in that it al- 
lows the measurement strengths to be fully tunable 
by adding more kicks, as well as the easy inclusion 
of unitary operations by simply waiting a fraction 
of a Rabi period. 

• Kicked QND measurements may be implemented 
in experiments with several advantages. The kicks 
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may be accomplished with a pulse generator on a 
QPC measuring a DD, and the waiting time be- 
tween kicks gives external decision circuitry the 
needed time to process the data in order to do 
real-time feedback. Also, the pump variation in- 
troduced in R.ef. I 111 removes the uninteresting back- 
ground signal of the measurement, and just gives 
the bare output signal as either positive or negative 
pumped current. 

The purpose of this paper is two-fold. The first topic 
is formal: To start with the well-known POVM approach 
to generalized measurements, and derive the quantum 
Bayesian formalism from it, starting with a scattering 
detector. The detector physics allows a natural transla- 
tion of the abstract POVM formalism into physical pro- 
cesses, and the quantum Bayesian formalism is recovered 
in the weak coupling limit. After discussing kicked-QND 
measurements, wc show how both kicked measurements 
and unitary operations may be recast in terms of con- 
formal maps, and demonstrate a close parallel with the 
mathematics of the special theory of relativity. 

The second topic is physical: the formal results are 
applied to make predictions using conditional operations 
with real-time feedback: (1) We derive an algorithm to 
deterministically produce an arbitrary pure state after 
a (random) weak measurement using feedback. (2) We 
investigate the purification process under measurement 
and generalize Jacobs' qubit feedback protocol to speed 
up purification with feedback^ 

The paper is organized as follows. In Sec. II, we de- 
rive the quantum Bayesian formalism from POVMs ap- 
plied to a mesoscopic scattering detector in the weak cou- 
pling limit. Kicked QND measurements are reviewed in 
Sec. Ill, in the context of the quantum Bayesian formal- 
ism. In Sec. IV we introduce a stereographic projection 
representation, and rewrite the measurement dynamics 
as a stochastic conformal mapping. A close analogy to 
special relativity is also discussed. Sec. V combines 
kicked measurements with unitary operations, and calcu- 
lates measurement statistics. Sec. VI introduces condi- 
tional phase shifts in order to deterministically produce 
the same quantum state after a measurement. In Sec. 
VII, we investigate the purification process of any ini- 
tially mixed density matrix under kicked measurement. 
Sec. VIII contains our conclusions. 



II. DERIVATION OF THE QUANTUM 
BAYESIAN FORMALISM FROM POVM 

The formalism used in this paper is called the quantum 
Bayesian approachJ^ because it may be considered as a 
generalization of classical Bayesian inference. An anal- 
ogous approach to quantum measurement that is bet- 
ter known in the quantum information community has 
been given the unfortunate name of positive operator- 
valued measure (POVM) formalism. -'^^■-'^^ In this section. 



the quantum Bayesian formalism for a solid state system 
is derived from POVMs. 

Consider a bipartite system composed of A and B, 
where the states of B are expressed in the orthonormal 
basis \Q)b, and the states of A are expressed in the or- 
thonormal basis \j) A- A unitary transformation that en- 
tangles the states in A with the states in B is given by 

|V')A|0)B^^MQ|^.)A|g)i3, (1) 

Q 

where we consider an initial state \'iP)a m A, described 
by a density operator p^, an initial state |0)b in B, and 
have introduced the measurement operators Mq that are 
indexed by the states in B, and operate in A. The nor- 
malization of the states gives the completeness relation, 
MqMq = 1. Now make a projective measurement 
on B alone, and find the result Q. Any measurement of 
this kind may be described as a POVM in A. The prob- 
abihty of finding the result Q, called P{Q), is given by 

P{Q) = Tr(p^M^MQ), (2) 

while the outcome of this measurement prepares a new 
density operator of A, conditioned on the outcome Q, 
and is given by 

Pa = i — — • (3) 

Tr(p^M^MQ) 

This defines a mapping p'^ = $(pa) from density oper- 
ators to density operators, known as a super-operator, 
which is not unitary in general. 

While the above results are standard generalizations of 
projective measurement on A, the abstract formulation 
obscures how to practically apply the POVM to a spe- 
cific physical system. We now consider such a system in 
the solid-state: the quantum point contact (QPC), mea- 
suring the state of a double quantum dot (DD), in order 
to see how the POVM translates into physical language. 
The Coulomb interaction between the DD and QPC al- 
ters the transport properties of the QPC, and can thus be 
used to detect which quantum dot the DD electron occu- 
pies. The QPC is described with the help of a scattering 
matrix Sj that depends on the physical state of the DD. 
Following Averin and Sukhorukov^ the unitary evolu- 
tion of the total state during the scattering process is 
comprised of the state of an individual electron (system 
B) incident from the left electrode, \in)B, and the state 
of the DD (system A), a\l)A + /3|2)a, evolving as 

\tn)Bia\l)A + mA)^ (4) 
ain\L)B + h\R)B)\l)A + /3(r2|i)s + t2\R) b)\2) a, 

where \L,R)b represent scattering states of the QPC that 
have either been reflected or transmitted, and the trans- 
mission and reflection amplitudes tj and rj are elements 
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of the scattering matrix Sj : 

The state of the "logical" qubit of the DD, is now entan- 
gled with the "ancilla" qubit of the left/right position of 
the QPC electron, and this comprises the bipartite sys- 
tem. Using the evolved state Q , we can now read off the 
measurement operators, Mg, of this unitary operation in 
the {\0)a, |1)a) basis, 

and easily verify that M|^Ml -f M^M^, = 1 from proba- 
bility conservation. Counting the electron in the collector 
of the QPC gives a random outcome, Q = 1 if the elec- 
tron is counted, or Q = if the electron is not counted, 
and makes a projective measurement on the B part of 
the Hilbert space. Equation jSJl gives the probability of 
counting the electron (or not), 

P{1) ^ PiiTi + P22T2, P{0)= P11R1+ P22R2, (7) 

where py are the elements of the DD density matrix in 
the |1, 2) basis, Tj = \tj\^, Rj = jr^p, and Tj + Rj = 1. 
The density matrix of the DD qubit may be updated, 
given the outcome of the measurement with Eq. (PJ. If 
Q = 1, so an electron is counted, then 

P'll = Tipii/P{l), P22 = l-Pll' 
P'l2 = {P2l)* =ht*2Pi2/P{l) 

= Pl2e*'^^p'iiP^2/PllP22, (8) 

where ^= Arg{tit2); while if Q = 0, so an electron is not 
counted, or equivalently, a hole is counted, then 

p'll = Ripn/P{0), p'22 = l-p'ii. 
P'i2 = (P21)* ^ nr*p,2/P{0) 

= P12 e'^-^ p'iiP'221 P11P22, (9) 

where x= Arg(rir2). 

The results (|8I9|) have a natural interpretation as a 
quantum Bayes formula: The diagonal density matrix el- 
ements are interpreted as classical probabilities, and are 
updated according to the classical Bayes formula, while 
the off-diagonal elements have a more exotic rule,^^ Note 
that if the initial DD qubit state is pure, it remains pure 
after the measurement. This is because while the entan- 
glement enlarged the effective Hilbert space which would 
lead to decoherence if the entangled information went 
undetected, the measurement of the QPC electron col- 
lapses the 2-particle state back down to a different pure 
DD state. 

It is instructive to contrast the POVM procedure with 
the well known "decoherence" approach to quantum mea- 
surement in this most simple case. The decoherence ap- 
proach corresponds to explicitly averaging the elements 



of the density matrix over all possible outcomes of the 
detector. In this case, the two possible outcomes of the 
measurement (I8I9|) are used to obtain 

(p'll) - P(0)p'ii(0) + P(l)p'ii(l) = pii (10) 

for the diagonal elements (p22 = 1 ^ Pii), and 

(p'12) =P(0)p'i2(0) + F(l)p'i2(l) = (<i<^+rir*)pi2 (11) 

for the off-diagonal elements (p2i = Pi2)^ in agreement 
with Averin and Sukhorukovi— iiS The new off-diagonal 
matrix elements are reduced because \t1t2 + t'i'^2 1 ^ 1; 
resulting in effective decoherence, while the diagonal ma- 
trix elements are preserved. The predictive advantage of 
the quantum Bayesian approach comes from not averag- 
ing over the measurement results, but rather conditioning 
the quantum density matrix on the result obtained in a 
particular physical realization. 

The above POVM analysis is not difficult to extend 
to M "ancilla" qubits, or QPC electrons. The basis \Q) 
is now spanned by M qubits, each being projected to ei- 
ther or 1. Rather than find the probability of obtaining 
a given sequence of Os and Is in the output, it happens 
that it is sufficient to find the probability of just obtaining 
the total charge n = X^iHi Qi^ given M attempts, where 
Qi — (0, 1). In other words, sequence does not matter, 
only the total number of counted electrons. This map- 
ping is illustrated in Fig. 1, where the ancilla outcome 1 
or 0, is mapped respectively into either counting an elec- 
tron, or not counting an electron in the current stream. 
The generalization of the quantum Bayesian rules H7I9|) 
for M ancilla qubits is done by replacing the success prob- 
abilities Tj , and the failure probabilities Rj = 1 — Tj by 
the probability P{m^ to measure m electrons in M 
attempts, under the condition that the qubit is in state 
\]), 

Pim, M\j) = (^^) Tp{l ~ T,)^^-™, (12) 

which is the binomial distribution. 

While this is sufficient for the diagonal matrix ele- 
ments, the off-diagonal matrix elements get a different 
phase shift after each electron is measured. The total 
phase shift is <I>(m, Af ) — Mx + ™(^ — x), which is de- 
terministic if ^ = X- The symmetric quantum point con- 
tact has the property that ^ — x = 0(Si4 so this sim- 
plification will be made in the rest of the paper. For 
many electrons M, the current is determined by m/M. 
The physically relevant weak-coupling (weakly respond- 
ing) limit corresponds to {Ti - T2) / {Ti + T2) < 1. Ap- 
pealing to the central limit theorem, we treat the detec- 
tor shot noise in the Gaussian approximation with little 
lost information. The QPC detector is efficient, in the 
sense that no information about the DD qubit is lost in 
the noisy current outputi^ii^ This is analogous to saying 
in the logical language that no ancilla qubit was left un- 
projected, and that the unitary operations did not hide 
any qubit information in the phase of the ancilla qubits 
that is destroyed after projection in the left/right basis. 
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FIG. 1: Weak entangling of a logical qubit with many ancilla qubits, followed by projective measurement on the ancilla qubits 
can be mapped onto the measurement of a quantum double dot by the transport electrons of a quantum point contact. Two 
given measurement realizations are shown, depending on whether the logical qubit state is either |1) or |0). The random ancilla 
result 1 is mapped onto measuring an electron in the current collector (denoted with a filled circle), while the random ancilla 
result is mapped onto not measuring an electron in the current collector, or equivalently, measuring a hole (denoted with an 
empty circle). The fact that the states |1) and |2) of the logical qubit alters the probability of projecting the ancillas to or 1, 
allows a weak POVM measurement on the logical qubit. 



One manifestation of an efficient detector in the de- 
phasing approach, is that the measurement rate coin- 
cides with the measurement-induced dephasing laXem^ 
Ref. related the measurement rate to one of the 
Renyi entropies, or the statistical overlap Oi^2 = 
X;,„[P(m,M|l)P(m,M|2)]i/2^ and showed that for the 
symmetric QPC it coincides with the reduction of the 
off-diagonal matrix elements. It is one feature of the 
quantum Bayesian approach that this particular measure 
comes out in a natural way. To see this, we first note a 
general property of the density matrix, \pi2\ < ^/PiiP22, 
that simply comes from the density matrix eigenvalues 
being bounded between zero and one. Next, we consider 
an initially pure state {\pi2 \ — ^/piTpii) and notice that 
the elements of the density matrix after a measurement, 
p', obey the relation 



Pl2 
Pl2 



< 



l p'iiP'22 
P11P22 



(13) 



The above relation is valid for every given measurement 
outcome, so it is also valid after averaging over the distri- 
bution of results, {0{m)) = X^m ^("^' -^^)'-^("^)' where 
O is any observable, and P{m,M) — piiP{in, M\l) + 
P22P{m, M\2). Using the classical Bayes rule for the di- 
agonal elements, taking O — \p'i2/ Pi2\, and the fact that 
\{^)\ 1^ we obtain the generalized efficiency rela- 

tion 



P12 
P12 



< 



y/P{m, M\l)P{m, M\2). (14) 



Notice this relation is quite general, as no particular up- 
date rule for the off-diagonal matrix element has been 
invoked. Therefore, a detector reaching the upper bound 
(|14|l can naturally be called ideal, or 100% efficient. (For 
the appropriate definiti on of e fffciency for an asymmetric 
detector see, e.g., Refs. ll5l2lL l 



Thus far, we have focused only on the dynamics of the 
measurement process, and have neglected the Hamilto- 
nian evolution of the DD qubit. This evolution rotates 
the quantum state, and continually changes the effective 
measurement basis, which typically ruins the desired con- 
tinuous measurement. The way to get around this both- 
ersome detail is with QND measurements, the subject of 
the next section. 



III. KICKED QND MEASUREMENTS 

The unifying theme behind all QND schemes is to cou- 
ple the measurement apparatus to the qubit with an op- 
erator that is an approximate constant of motion of the 
measured quantum systemi^ In this way, the detector 
only measures the state in the desired fixed basis, and the 
internal quantum dynamics that would otherwise spoil 
the desired measurement is circumvented. The specific 
scheme we employ in this paper is that of kicked QND 
measurements, introduced by V. Braginsky et al7 and 
K. Thorne et alMi for the harmonic oscillator. In Ref. 9), 
this idea is introduced for two-state systems by making 
an analogy to a cat playing with a string that moves 
in a circle. In the kicked QND mode, the cat sits in 
one spot waiting for the string to come to it, and only 
then bats at it. The motion in a circle comes from the 
simple Hamiltonian evolution of a two-state system. If 
H — 60-2/2 + Aax/2 is the qubit Hamiltonian, where A 
is the tunnel coupling energy, and e is the energy asym- 
metry, then unitary evolution for a time t is given by 

U = exp[-it{e(T, + A(Tx)/2] = lcos(M/2) (15) 
- ia^{elE) sin(M/2) - ia^iA/E) sm{Et/2), 



where E = Ve^ -|- A^, and h — 1 throughout the paper. 
From the perspective of the qubit, the measurement ap- 
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FIG. 2: (color online). (After Ref. ,11). Visualization of the 
kicked QND measurement scheme. A voltage pulse is applied 
to the quantum point contact on a time scale tv <C t,, fol- 
lowed by a quiet period of zero voltage bias, lasting for a 
Rabi oscillation period t,, followed by another pulse, and so 
on. The up/down variation is depicted, where the kicks come 
every half period, and the sign of the voltage pulse alternates 
with every kick. In this scheme, qubit read-out is by sim- 
ply measuring the sign of the current, and corresponds to an 
elementary quantum pump. 



paratus only measures at approximately discrete points 
in time. In this reduced problem, by choosing the wait- 
ing time between kicks to be = 2'k/E (or some integer 
multiple n thereof), the unitary evolution IjlSI) becomes 
U (—1)". The operator we want to measure is then 
static in time, and is thus a QND measurement. (The 
evolution is also simple if n is a half-integer, especially if 
e = 0). However, from the point of view of the detector, 
the on/off pulse lasts much longer than any detector time 
scale, so many electrons pass through the QPC. If tq is 
the time scale of the QPC electron correlation, ry the 
time scale of the pulse duration, and Tq is the Rabi oscil- 
lation period, then the considered time scale ordering is 

To < TV < Tq^ 

A pump variation on kicked QND measurements was 
given by Biittiker and the authors in Ref. Illl where in- 
stead of giving the same kick every Rabi oscillation, the 
experimentalist gives a sequence of voltage kicks to the 
QPC with a pulse generator, alternating in sign, every 
half oscillation period (see Fig. 2). In this scenario, 
we have shown that if e = 0, qubit readout is accom- 
plished by pumping current: the kicks provide one AC 
current source, and the dynamics of the qubit provides 
another (intrinsically quantum mechanical) AC current 
source, that nevertheless causes a net DC current flow in 
the QPC j^^i^"^ There are two limiting cases the system is 
driven into: Either the qubit oscillations are of the same 
phase with the pump oscillations, pumping positive cur- 
rent, or the qubit oscillations arc of opposite phase with 
the qubit oscillations, pumping negative current. 

To characterize the result of each measurement kick, 
the parameters of the measurement process with an 
ideal QPC detector are specified by the currents, h 
and I2, produced by the detector when the qubit is 
in state |1) or |2), and the detector shot noise power 
Si = e/(l - T) (where T is the transparency).^ The 
typical integration time needed to distinguish the qubit 



signal from the background noise is the measurement 
time Tm — ^Si/{Ii — 12)^- Shifted, dimensionless vari- 
ables may be introduced by defining the current origin 
at /o = [h + /2)/2, and scahng the current per pulse as 
I — Iq ~ x{Ii — /2)/2, so /i_2 are mapped onto x — ±1. 
The weak static coupling (per pulse) between QPC and 
DD implies that the kick duration ry is less than the 
measurement time Tm- We take x to be normally dis- 
tributed with variance D — T^i/tv- The typical number 
of kicks needed to distinguish the two states is D, where 
we assume D ^ 1. 

The measurement result T after N kicks is 



1 



N 



(16) 



and we seek the conditional probability distribution 
P{X,N\p) of measuring the dimensionless current X, 
starting with a given density operator p prepared be- 
fore the first kick. The functions P(I, N\j) are defined 
as classical probability distributions of the current with 
mean I = ±1 (if j = ±1) and variance cr"^ = D/N, the 
Gaussian equivalent of H12|l . 



P(X,iV|l) = 
P(X,iV|2) = 



1 



^2ttD/N 
1 



exp 



exp 



(I -If 



2D IN 
{1+lf 



2D/N 



(17) 



and the notation Pj{xn) is adopted for the j — 1,2 dis- 
tributions of the nth kick. The probability density of 
measuring the result Xn after one kick is determined by 
the state of the qubit just before the measurement, and 
is given by the analog of {Tj), 



P(a;„) = pii'Pi(x„) + pt^P2{xn)- 



(18) 



The density matrix of the qubit is updated based on the 
information obtained from the measurement that just oc- 
curred. This is done with the quantum Bayesian update 
rules, that defines a non-unitary quantum map directly 
analogous to Eqs. (|8I9|I . 



(n+l) 
Pll 



p^'^ PliXn) + P^22 P^iXn)' 



^ Pll 



Jn+1) _ 
P12 — 



(19) 



This quantum map is a probabilistic, non-unitary relative 
of the unitary maps studied in kicked quantum chaosi^ii^ 
The advantage of QND measurement in the Bayesian ap- 
proach is seen by using Eqs. H19|) to express the condi- 
tional probability density P(a;„) in terms of the result of 
the preceding kick Xn-i- It follows from (|18I19|1 that 

lr'^^l(a;»-l)Pl(a;„) + p'2r'^P2{Xn-l)P2{Xn) 



P{Xn) 



Pi 



P{Xn-l) 



(20) 
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This recursive relation helps in the calculation of the (un- 
conditional) probability distribution P{T, N\p) of finding 
current T, starting with the density matrix p, after N 
kicks, given by 

P(I, N\p) = / n dXnP{Xn) ^^/^ ) ■ (21) 
n=l \ i=l / 
I 



Each application of (|20|l generates a denominator that 
cancels the probability density immediately preceding it 
in inj. Making N iterations of gives 



P{I,N\p) 



N 

I W dxn [pii-Pi(a;i) . . . Pi (a; at) + P22P2{xi) ■ 



N 



. P2ixN)] 5{I~Y1 ^^/^) = -^ll^^^' N\^)+P22P{I. N\2), 

(22) 



4=1 



where pii,p22 are the diagonal matrix elements of the 
original density matrix, and the N Gaussians compose 
to form one Gaussian with a variance N times smaller. 
As N is increased, the two qubit states can be distin- 
guished with greater statistical confidence, and eventu- 
ally the distributions limit to delta-functions, giving ei- 
ther I — 1 with probability pn , orX — —1 with probabil- 
ity P22- A one-sigma confidence is obtained when N = D, 
as previously stated. 

It is worthwhile to point out several features of the 
above QND measurement. First, N weak measurements 
simply compose to make an iV-times stronger measure- 
ment. Second, the QND measurement output only in- 
volves the diagonal density matrix elements. It is for 
this reason that the output of a quantum nondemoli- 
tion measurement is equivalent to noisy classical mea- 
surement, where the detected "classical probabilities" are 
given by the diagonal density matrix elements in the 
preferred measurement basis. In spite of the classical 
nature of the detector output, the qubit state prepared 
after the N measurements can be deduced from the out- 
come of the random variable X. To characterize the post- 
measurement density matrix, we note another recursion 
relation from IjlQII . 

pir^Vp^'^ = [p[l^/p'2^^][PlM/P2{Xr.)] 

= [pliVp^2^]exp(2a:„/i?). (23) 

This result may be composed N times, and the update 
of the off-diagonal matrix element H19|l follows from (|23|l . 
Using the definition X = {^/N)^^^^Xi, the measure- 
ment of a given random current I prepares a new density 
matrix of the DD, 
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P 



pii el + P22e 1 



Pii e 

Pl2 



7 



Pl2 
P22 e~ 



(24) 



where 7 — XN/D is the rescaled measurement result, 
named the rapidity of the measurement, for reasons given 
in the Sec. IIVI The conditional quantum dynamics of 
Eq. (|24|l is illustrated in Fig. 3, for all pure states, and the 
density matrix is parameterized as p = (\ -\-^^Xiai) /2, 
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FIG. 3: (color online). (After Ref. Q. The conditional 
evolution of all initial pure states, represented on the Bloch 
sphere, under continuous QND measurement by an efficient 
detector. From (a-e), the rapidity of the measurement is 
7 = NX/D — ( — 1, — .5, 0, .5, 1) respectively. As the detec- 
tor obtains more information about the quantum state, we 
can with greater statistical certainty distinguish the post- 
measurement quantum state, so the Bloch sphere is more and 
more red (|1)) or blue (|2)), depending on the value of the 
rapidity measured. The conditional evolution of several rep- 
resentative states is also indicated with black arrows. The 
view-point is parallel to the equator of the Bloch sphere, so 
|1) North pole, and |2) -> South pole. 

so {X, Y, Z) give coordinates on the Bloch sphere. The 
X and Y behavior follows from Z, which is in turn condi- 
tioned on the detector output X, so the sphere is colored 
according to the conditional evolution of Z, 



Z' 



sinh "/ + Z cosh 7 
cosh 7 + Z sinh 7 



(25) 



If the rapidity 7 is positive, then states are "attracted" 
toward the North pole, while if the rapidity 7 is nega- 
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tive, then states are "attracted" toward the South pole. 
As the rapidity grows increasingly positive or negative, 
we become more confident which state the qubit has con- 
tinuously collapsed to, but this also depends on the initial 
state. The conditional evolution of several representative 
states is indicated with black arrows. 

While the whole point of the kicked QND proposal was 
to effectively turn off the qubit unitary evolution while 
the continuous measurement is taking place, a much 
more interesting situation arises when continuous (non- 
unitary) measurements are combined with controlled uni- 
tary rotations. The kicked measurement set-up provides 
a simple way of generating a single-qubit rotation: wait- 
ing. Rather than spacing the pulses by a full Rabi os- 
cillation as previously described, we choose to wait some 
fraction r of a Rabi oscillation, twait = fTq, that defines 
a phase shift (j) — 2'Kr. A single-qubit unitary operation 
(expressed in the z-eigenbasis). 



U = 



a -b* 
b a* 



(26) 



may be executed by choosing r, such that a — cos((/)/2) — 
iie/E) sin(0/2), and b = -i{A/E) sm{(j)/2), by using the 
Hamiltonian evolution of the qubit, Eq. H15|l . Varying r, 
any point may be reached on a circle on the Bloch sphere, 
which is fixed by e and A. In order to reach any pure 
state (up to an overall phase) by Hamiltonian evolution 
starting with a pure state, the qubit asymmetry e should 
also be varied between the kicks with gate voltages. 



IV. POVM MEASUREMENT AS A 
STOCHASTIC CONFORMAL MAP 

Before considering specific examples of weak measure- 
ment, combined with unitary operations, we first refor- 
mulate the above set of weak measurements (|18ll9ll , and 
unitary operations H26I) in an important special case: 
where the initial state is pure, the detector is efficient 
(as considered in this paper), so the post-measurement 
state is also pure. The initial arbitrary DD state is 
defined as = a\l) + /3|2), with density matrix ele- 
ments pii = |ap, pi2 = al3*, P21 = PI2, P22 = = 
1 — pii. Represented as coordinates on the Bloch sphere, 
{X, Y, Z) , both measurements and unitary operations 
leave the state on the surface, X^ + Y"^ + Z"^ = 1. 

Now make a stereographic projection of the Bloch 
sphere onto the complex plane, with the complex vari- 
able (, defined as 



C = P12/P22 = {X + iY)/{l - Z) = a/p. 



(27) 



The South pole of the Bloch sphere is identified as the 
origin of the ( plane, while the North pole is identified as 
00 on the C plane. Translating the unitary operation H26() 
on the qubit into an operation on the complex variable 
C, we find the conformal mapping. 



aC~b* 
bC + a*' 



(28) 



known as a Mobius transformation. The group algebra 
of unitary rotations maps onto the group algebra of con- 
formal Mobius transformations. Reference i29l points out 
that this property can be used to demonstrate operator 
product identities for one qubit. 

Translating the non-unitary Bayesian update equa- 
tions for the density matrix H19|) . as an operation on the 
complex coordinate ^, we find the following stochastic 
conformal mapping: 



C = C ^Pi{x)/P2{x) = C exp(x/i?). 



(29) 



This conformal mapping is simply a random scale trans- 
formation with two fixed points: one at (the South pole, 
state 1 2)) and the other at 00 (the North pole, state |1)). 
The random variable x is chosen from the probability 
distribution (|18|) . which translates to 



P{x) 



cPi{x) + icr'P2{x) 



(30) 



Thus, any sequence of weak measurements, combined 
with unitary operations can be translated into repeated 
conformal mapping. We note that after an arbitrary se- 
quence of weak measurements and unitary operations, 
the definition C = a/ (3, together with the normalization 
of the state, |ap-|-|/3p = 1, immediately allows the wave- 
function to be read off (up to an overall phase). The 
inclusion of asymmetric measurements, where the phase 
shift in Eq. (O is kept for a broader class of scattering 
matrices may also be easily included. This phase shift 
has the effect of twisting the Bloch sphere proportionally 
to the value of current measured, which is equivalent to 
including phases in the scale factor. 



C' = C expiiOi) exp[a;(l + ie2)/D], 



(31) 



where 9i , 62 correspond to the continuous limit of the 
total acquired phase shift $(m, M) = Mx + — x), 
described in Sec. |n] The mapping (|31() is obviously still 
conformal. 

If we momentarily let a; be a deterministic variable, 
then the set of Mobius mappings form a group, and the 
set of deterministic scale transformations form a differ- 
ent group. A natural question that arises is what (if 
any) group is described by the composition of H28|l , with 
(|29|l ? Amusingly, the answer is provided by the special 
theory of relativity. Consider a relativistic four-vector 
(AT, y, Z, T) . As is well known, any element of the Lorentz 
group may be produced by making a spatial rotation, 
followed by a boost in the (say) Z direction, followed by 
another spatial rotation. The boost from (A, F, Z, T) to 
(A', F', Z' , T') in the Z direction may be described as a 
hyperbolic rotation 



( cosh 7 sinh 7 
lsinh7 cosh 7 

(1/2) log 



1 



X' = A, Y' = Y, 



(32) 
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where the rapidity 7 is introduced in terms of the velocity 
parameter v, the physical velocity measured in units of 
the speed of light, c—1. 

To connect this to spinor formalism, we follow the dis- 
cussion in Penrose and Rindlcr,— and define a Hermitian 
coordinate operator, 

^_ f T + Z X + iY\ , . 

^-[x~iY T-Z )■ ^-^-^^ 

Translating the boost H32() into an operation on the co- 
ordinate operator C yields 

C'^ACAT, A^(^^ (34) 

By fixing T = 1, the celestial sphere X'^ + Y"^ + Z"^ = \ 
is defined. The celestial sphere is then stereographi- 
cally projected, defining the complex variable C = (-^ + 
iY)/{l — Z). In the complex plane, the boost is simply 
a scale transformation, 

C' = Cexp7. (35) 

To extend the analysis to N boosts with rapidities 7^, 
the mapping simply composes the N scalings, to pro- 
duce a boost with rapidity 7 = Y^^=i li- This conformal 
mapping is identical with (|29|l . the analogous quantum 
measurement composition, if 7 = NT / the quantum 
measurement parameter, is identified with the rapidity 
of the boost. Furthermore, any spatial rotation of the 
sphere may be interpreted as a unitary operation on the 
Bloch sphere, which projects to the Mobius mapping (up 
to an overall phase). Thus, the group described by the 
composition of H28|) . with (|29|) is the (restricted) Lorentz 
group. 

The difference with the relativity analogy comes when 
we recall that 7 = X^i^*/^ ^ random variable. The 
distribution of this random variable l|3U|) . explicitly de- 
pends on the "space-time" coordinates C, and thus breaks 
the Lorentz invariance by introducing a preferred refer- 
ence frame, the Z-axis, with Z = ±1 as the attracting 
fixed points. From the quantum measurement point of 
view, this is a consequence of choosing to measure along 

I 



and the natural Hamiltonian dynamics performs the 
unitary operation (|15I26|I . One interesting feature of 
the result H38I39|1 is that the outcome of the first A^i 



the Z-axis. Therefore, for pure states, the mapping (|29l 
I30|l may be viewed as a stochastic Lorentz semi-group. 



V. COMBINED WEAK MEASUREMENTS AND 
UNITARY OPERATIONS 

After having separately described weak measurement 
and unitary operations, we now combine them. Consider 
an experiment, where A'^i kicks are made, followed by a 
single qubit unitary operation U (produced by inserting 
a dislocation into the pulse sequence), followed by N2 
kicks. The measurement results Xi and T2 are defined as 

^ ^ N1+N2 

Xi-— ^x., J2 = _ ^ X,. (36) 
^ j=l ^ 1=^1+1 

We seek the normalized probability distribution 
P{Ti, Ni;T2, N2) of finding current Ti after iVi kicks, 
and X2 after N2 subsequent kicks. This distribution may 
also be interpreted as joint counting statistics. 

The analysis from Sec. IIIII indicates that after the first 
Ni kicks, the measured current Xi will occur with a prob- 
ability given by (|22|) . and prepares a post-measurement 
density matrix p', given by (|24|l . described with the ra- 
pidity of the measurement j = XiNi/ D. The subsequent 
unitary operation U H26|l . (characterized by a phase (p) 
rotates the post-measurement density matrix, 

p"°" = Up'U^ (37) 

The following set of N2 kicks start with the density ma- 
trix l|37|) . and continue to measure in the z-basis as be- 
fore. Equation H22I) may be applied again with the mod- 
ified initial density matrix (|37|l to deduce the (uncondi- 
tional) probability distribution of finding result Xi after 
A^i kicks, and result T2 after N2 kicks, 

P{Ji,Ni;J2,N2) - P(Xi,7Vi|p) X P(X2, A^2|p"°"), (38) 

where P{Ti, Ni\p) is given in 1221), and the new density 
matrix elements are given in terms of the phase (j) and 
rapidity 7 as (we let e = for simplicity), 



(39) 

I 

measurements, Ti, appears in the expression involving 
the variables of the second set of kicks. This immedi- 
ately implies that the distribution does not factorize. 



cos2(0/2)piieT -hsin2(0/2)p22e"'^ - 2 sin(0/2) cos(0/2) Impi 
PiieT + p22e^^ 



Pll — I P22 — ^ Pll 1 



now _ Repi2 + (z/2)sin(/) (pne^ - p22e ^) + i cos Impiz ^ow _ / ncwN* 

Pi2 — r~~ zr^ J P21 — \Pi2 ) 

piie^ + p22e T 
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P(Ji,7Vi;X2,A^2) Pi{Ii,Ni)P2{l2,N2). The effect 
comes from the first set of measurements preparing a 
given density matrix of the DD, which affects the results 
of the next set of measurements. From the distribution 
(|38|1 . the average current in each interval, as well as the 
correlation between the two may be calculated: 

(li) = Pii - P22, 

{I2) = (pii -p22)cos0-2exp(-A^i/2D)sin0Im/9i2, 
(XiX2> = COS0. (40) 



Also, as A^i, N2 are taken to infinity in (|38I39|I . the distri- 
bution Ppm(Ii,I2) from making simple projective mea- 
surements on the DD is recovered, 



P 



PM 



= pnCos\cf>/2)6{Ii-l)S{l2-l) 



Piisin2(0/2) S{Ii 
P22 sin2(0/2) S{Ji 

P22 0082(0/2) S{Ii 



l)S{l2 + 1) 

1)<5(I2 - 1) 
-1)<5(X2 + 1). 



(41) 



It is now straightforward to generalize the result H38I39|) 
to any number of m—1 dislocations in the pulse sequence, 
each of which has a phase shift of 4>k (and now e is arbi- 
trary), 

m 

P{{I„N,}) = n [p['l^P{Ik,N,\l)+pg^P{I,,N,\2) , 
fe=i 

(42) 

and each density matrix pC^+i) is defined in terms of the 
density matrix p^''^ after the previous dislocation. 



where = p^ii e.'^'' + p^22^~^'' ^ rapidities are 7^ = 
X^Nk/D, the initial density matrix before the first kick 
is p*-^^ , and the matrix U^: has elements 

(Ufe)ii = cos(0fe/2)-z(e/S)sin(0fe/2), 

(Ufe)i2 = -z(A/£;)sin(0fe/2), 

(Ufe)2i - -i(A/£;)sin(0fe/2), 

(Ufe)22 = cos(0fc/2) + i(e/S)sin((/.fe/2). (44) 

Reference applied these results to violate a Bell in- 
equality in time. 



u 



(43) 



VI. CONDITIONAL PHASE SHIFTS AND 
FEEDBACK PROTOCOLS 

In the preceding Section, the phase shift (/) was cho- 
sen beforehand, independently of the result Xx. We can 
now use the information gained in the first A^i measure- 
ments, and make a conditional phase shift, pending the 
outcome of the random variable X\ . This is essentially a 
feedback protocol that the experimentalist can choose to 



execute, defined by a function (j){Xi), so a different phase 
shift is assigned to every possible random outcome of the 
continuous measured current. 

Kicked QND measurement provides a realistic mecha- 
nism for implementing general qubit feedback protocols. 
The reason for this is two-fold: First, as seen in the previ- 
ous section, any combination of weak measurements and 
unitary operations may be accomplished with a sequence 
of voltage pulses to the detector. Second, the feedback 
circuitry must take the result obtained from the measure- 
ment, execute logical operations, and command the ex- 
perimental apparatus to do something it otherwise would 
not have done (like make a given phase shift). The intrin- 
sic waiting time between the kicks provides the needed 
time delay for all of the above to take place. 

We now explicitly find the feedback protocol 4i{Xi) to 
take a given pure state to any desired pure state after 
a weak measurement. As the simplest case, consider 
any pure state on the Z ~ Y great circle of the Bloch 
sphere {X = Q,Y^ + = 1), and a symmetric qubit, 
6 = 0. Both Hamiltonian evolution, and weak measure- 
ments (according to (|24|l ) do not take these states out of 
the Z — Y great circle. Therefore, knowing the outcome 
Xi , a conditional phase shift may be applied to determin- 
istically prepare any quantum state on the Z — Y great 
circle of the Bloch sphere. For dcfiniteness, we choose 
to shift to the state |1). This choice has the advan- 
tage that after the N2 ^ D measurements, the current 
will be X2 = 1 deterministically. The parametrization 
{Y,Z) — (— sin 6*, cos 0) of the initial state is chosen, so 
that if no measurement is made, the shift to the North 
pole may be done with a phase shift (j) = 9. The result 
is applied by setting p^™ = 1, and solving for (p as 
a function of the rapidity 7, to find 



tan((/)/2) =tan(6l/2) exp(-7). 



(45) 



This answer interpolates between two extreme 
strategies!^ (1) If no measurement is made, just 
make the desired phase shift, (j) = 9. (2) If a projective 
measurement is made, either do nothing if Xi = 1, or fiip 
the state by applying the phase shift = tt if Xi = — 1. 
The asymptotic limits in the later case may be obtained 
by expanding the inverse tangent to obtain 



0/2 



tan (0/2) exp(-7) 
,f «ign(0)-t^ 



if 7 > 1, and 9 ^ ±n, 
if 7 < -l,and 6* 7^ 0. 

. 

The above real-time feedback proposal is experimen- 
tally promising in the kicked scheme. However, it de- 
mands fast time resolution and feedback circuitry. An 
experimentally simpler proposal to verify the above pro- 
tocol is to make many realizations of weak measurement, 
phase shift, weak measurement, where the phase shift is 
chosen randomly in each realization. After the run is 
finished, the data record may be reviewed, and all in- 
stances of phase shifts where condition H45|) is approx- 
imately satisfied are post-selected. In this data subset, 
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the prediction is that the following set oi N2 ^ D kicks 
will deterministically find X2 = 1. 



VII. PURIFICATION OF INITIALLY MIXED 
DENSITY MATRICES BY WEAK 
MEASUREMENT 

Under repeated weak QND measurements, eventually 
all states collapse to either |1) or |2), including mixed 
initial states. However, the states |1) and |2) are both 
pure, and therefore if the initial state is mixed, a purifica- 
tion occurs during the measurement process^ This phe- 
nomenon is especially counter-intuitive from the point 
of view of the dephasing approach to quantum measure- 
ment. Jacobs has shown that the average purification 
in a given time can be increased by the use of continu- 
ous feedback.^** Jacobs' protocol is somewhat counterin- 
tuitive for the qubit: always use Hamiltonian evolution 
to rotate the state to Z = 0, i.e. perpendicular to the 
measurement axis. The purpose of this section is (1) To 
show how this idea can be easily implemented for our 
set-up, (2) To demonstrate that the "equatorial plane" 
protocol {i.e. Z = 0) is also optimal for kicked QND 
measurements which have a continuous output of tunable 
measurement strength, and (3) To generalize Jacobs' no- 
feedback purification solution to any initial density ma- 
trix. 

It is well known that any unitary operation preserves 
the purity (or entropy) of the state. It is interesting to 
note from (|19|l . that during weak measurement there is a 
different preserved physical quantity, that we name the 
murity. For the qubit, the purity V and the murity M 
are defined as 



V = 



M ={X^ + Y^)/{1- Z^). (47) 



If the purity V — 1, then the murity M — 1, reflecting 
the statement made in Sec. m that if the initial state 
is pure, the post-measurement state is also pure. We 
note that V may be expressed in terms of M and Z by 
V M{1 - Z'^) + Z^. After one kick, the change in the 
purity, AV, (or purification) is given by 



AV = V' - r = {1- M)[{Z')^ - Z\ 



(48) 



where we have used the fact that murity does not 
change during measurement. Application of the quan- 
tum Bayesian update rules l(T^ yields Z' = [piiPi(a::) — 
P22^2(a:)]/(piiPi(x) + P22P2{x)) (also given in 
where x is the measurement result, so the purification 
is 



AV ^{l-V)[l 



[cosh{x/D) +ZsiTi\i{x/D)Y 

'(49) 

Several observations are in order: First, li V = 1, the 
purification AV is automatically 0, while the first (de- 
terministic) factor is maximal liV = Q. Second, if x = 0, 



the second (random) factor is zero, so there is no pu- 
rification, which corresponds to no gained information. 
Finally, the first factor is between [0, 1], while the second 
factor could be negative or positive, implying that either 
purification or further mixing is possible in a given run. 

The average purification is given by averaging 1)49(1 over 
the distribution of x, given in Eq. IIHI), to yield 



{AV) = {l-V)[l- f{D,Z)l 



(50) 



where 



I{D.Z) 



dx 



exp(-a;V2-D) 



, ^J2ttD cosh{x/D) + Z sinh(a;/i:') ' 

(51) 

It is straightforward to check that < / < 1, so 
there is nonnegative average purification for all den- 
sity matrices. ^'^ Changing variables to 7 = x/D, it is 
also straightforward to check asymptotic limits. Tak- 
ing D ^ 0, the projective hmit, /(O, Z) = is recov- 
ered, so {AV) = 1 — V, implying that the final state is 
pure with unit probability. The opposite limit, D —^ 00, 
corresponds to an vanishingly weak measurement, so 
/(oo, Z) = 1, or {SV) ~ 0, giving no purification. 

The results (|50I51|I allow us to find the optimum aver- 
age purification strategy for one kick. The best strategy 
on average is to rotate the density matrix to where the 
purification (|50|l is maximum. These point (s) may be 
found by maximizing {AV) on the Bloch ball, under the 
constraint that D and V are fixed. The V constraint sim- 
ply reflects the fact that unitary operations do not alter 
the purity. This problem is equivalent to minimizing / 
with respect to Z, by solving df /dZ = 0, which leads to 
the equation 



c?7 



exp(— i?7^/2) sinh7 
(cosh 7 + Zs'mh'y)^ 



= 0. 



(52) 



The solution is immediate because at the point Z — 0, 
the integrand is odd, so the integral is zero. The fact 
that / is minimized at Z = is seen by noting that 
the integrand of d^f/dZ^{Z = 0) is nonnegative. There- 
fore, the average purification is maximized by applying a 
phase shift after the first measurement that rotates the 
qubit to the equatorial plane of the Bloch ball before the 
measurement. This result is in agreement with Jacobs, 
who considered purification from a two-outcome POVM 
of variable strength (and also the stochastic Schrodinger 
equation limit) li^ Our approach is from the complimen- 
tary perspective of a continuous outcome measurement 
of variable strength. 

Results (|50l51ll have a simpler form in the large D 
limit, for very weak measurements. The average purifi- 
cation {AV) and the noise in the purification, {{AV}'^), 
are given to leading order in D^^ as 



{AV) = {1-V){1-Z^)/D, 
((AP)2) ^i{l-V)^zyD. 



(53) 



In order to compare purification with and without feed- 
back, we first consider Jacobs' feedback protocol Z = at 
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every time step. Equation (|53() implies that for this feed- 
back protocol the purification noise vanishes. Therefore, 
the dynamical purification is described by the determin- 
istic rate equation dV/dN = (1 — V)/D. Solving the 
equation with initial condition Vq yields^"' 

{Vn) - 1 + (Po - 1) eM-N/D), (54) 

showing an exponential approach to a pure state, with 

I 



rate D 



In the no feedback case, the purity after N kicks may 
be found from the murity relation V' = M[l — [Z'Y] + 
[Z'Y' , where Z' is given in (|25|l . Averaging this relation 
over the distribution H22|l yields the average purity after 
N kicks,^ 



D 
2ttN 



exp 



N 
2D 



—— / d7exp 



Dj'^\ M(cosh7 + Zsinh7)2 - (Af - l)(sinh7 + Zcosh7)^ 



2N 



cosh J + Z sinh 7 



(55) 



r 



After some manipulation, the above integral expression 
may be simplified to 



(Vn) = 1 - {M - 1){Z^ - 1) 



D 
2nN 



exp{-N/2D) 



d'y 



exp{-D-f^/2N) 
cosh 7 + Z sinh 7 



(56) 



For large N/D, the dominant dependence comes from 
the term outside the integral, and the N/D dependence 
inside the integrand may be neglected. In this limit, the 
purity may be approximated as 



(Vn) ~l-^T{l-M)^/^~^^/DJ(2^exp{~N/2D), 

(57) 

yielding an approach to purity with rate {2D) ^, half as 
fast as the feedback case, in agreement with Jacobs. The 
result 156|) generalizes Jacobs' no-feedback purification 
result to arbitrary initial states. Before concluding, we 
point out that Wiseman and Ralph have recently shown 
that the advantage of feedback for purification depends 
on how the question is formulatedi^^ If instead of asking 
about the average purification for a fixed time, we ask 
about the average time taken to reach a given purity, 
then the no feedback case is actually better. 



VIII. CONCLUSIONS 

The quantum Bayesian approach to the problem of 
quantum measurement has been derived from POVM 
formalism, applied to a mesoscopic scattering detector. 
By considering an elementary scattering event, measure- 
ment operators associated with the successful or failed 
detection of the electron in the current collector can be 
identified. We recover the quantum Bayesian formalism 
in the continuous current approximation. 

Kicked QND measurements have been analyzed within 
the quantum Bayesian formalism. We derive a quantum 
map representation that, while discrete in the time index. 



describes a sequence of weak measurements. Unitary op- 
erations (easily implemented by waiting a fraction of a 
Rabi period), together with kicked measurements, can be 
represented as a sequence of conformal mappings, where 
the unitary maps are deterministic, and the kicked mea- 
surement maps are stochastic. A close analogy exists 
between these quantum maps and the Lorentz transfor- 
mations of special relativity. 

We have calculated the measurement statistics asso- 
ciated with combined weak measurements, and unitary 
operations. These results are applied to find the feedback 
protocol that deterministically takes a given pure state to 
any other desired pure state after a weak measurement, 
using conditional phase shifts. 

Next, we have investigated the process of purification 
of mixed density matrices under kicked QND measure- 
ment. The concept of "murity" (the physical quantity 
that is preserved under measurement) has been intro- 
duced, and applied to calculate the change in the state's 
purity associated with a measurement. Purification with 
and without feedback has also been investigated. 

We stress that kicked QND measurements provides an 
experimentally viable way of implementing ideas in quan- 
tum feedback: Any combination of weak measurements 
and unitary operations can be accomplished by apply- 
ing a sequence of voltage pulses to the detector, and the 
intrinsic quiet time between kicks allows the necessary 
processing time for feedback to occur. 
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